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ABSTRACT 

We investigate the dependence of galaxy clustering on luminosity and spectral type 
using the 2dF Galaxy Redshift Survey (2dFGRS). Spectral types are assigned using 
the principal component analysis of Madgwick et al. We divide the sample into two 
broad spectral classes: galaxies with strong emission lines (‘late-types’), and more 
quiescent galaxies (‘early-types’). We measure the clustering in real space, free from 
any distortion of the clustering pattern due to peculiar velocities, for a series of volume- 
limited samples. The projected correlation functions of both spectral types are well 
described by a power law for transverse separations in the range 2<{a/h~^ Mpc)<15, 
with a marginally steeper slope for early-types than late-types. Both early and late 
types have approximately the same dependence of clustering strength on luminosity, 
with the clustering amplitude increasing by a factor of ~2.5 between L* and 4L*. At all 
luminosities, however, the correlation function amplitude for the early-types is ~50% 
higher than that of the late-types. These results support the view that luminosity, 
and not type, is the dominant factor in determining how the clustering strength of the 
whole galaxy population varies with luminosity. 


Key words: methods: statistical - methods: numerical - large-scale structure of 
Universe - galaxies: formation 


1 INTRODUCTION 

One of the major goals of large redshift surveys like the 2 de¬ 
gree Field Galaxy Redshift Survey (2dFGRS) is to make an 
accurate measurement of the spatial distribution of galaxies. 
The unprecedented size of the 2dFGRS makes it possible 
to quantify how the clustering signal depends on intrinsic 
galaxy properties, such as luminosity or star formation rate. 

The motivation behind such a program is to character¬ 
ize the galaxy population and to provide constraints upon 
theoretical models of structure formation. In the current 
paradigm, galaxies form in dark matter haloes that are built 


up in a hierarchical way through mergers or by the accre¬ 
tion of smaller objects. The clustering pattern of galaxies 
is therefore determined by two processes: the spatial dis¬ 
tribution of dark matter haloes and the manner in which 
dark matter haloes are populated by galaxies (Benson et al. 
2000b; Peacock & Smith 2000; Seljak 2000; Berlind & Wein¬ 
berg 2002). The evolution of clumping in the dark matter 
has been studied extensively using N-body simulations of 
the growth of density fluctuations via gravitational instabil¬ 
ity (e.g. Jenkins et al. 1998; 2001). With the development of 
powerful theoretical tools that can follow the formation and 
evolution of galaxies in the hierarchical scenario, the issue of 
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how galaxies are apportioned amongst dark matter haloes 
can be addressed, and detailed predictions of the clustering 
of galaxies are now possible (Kauffmann, Nusser & Stein- 
metz 1997; Kauffmann et al. 1999; Benson et al. 2000a,b; 
Somerville et al. 2001). 

The hrst attempt to quantify the difference between 
the clustering of early and late-type galaxies was made us¬ 
ing a shallow angular survey, the Uppsala catalogue, with 
morphological types assigned from visual examination of the 
photographic plates (Davis & Geller 1976). Elliptical galax¬ 
ies were found to have a higher amplitude angular corre¬ 
lation function than spiral galaxies. In addition, the slope 
of the correlation function of ellipticals was steeper than 
that of spiral galaxies at small angular separations. More 
recently, the comparison of clustering for different types has 
been extended to three dimensions using redshift surveys. 
Again, similar conclusions have been reached in these stud¬ 
ies, namely that ellipticals have a stronger clustering ampli¬ 
tude than spirals (Lahav & Saslaw 1992; Santiago & Strauss 
1992; lovino et al. 1993; Hermit et al. 1996; Loveday et al. 
1995; Guzzo et al. 1997; Willmer et al. 1998). 

The subjective process of visual classification can now 
be superseded by objective, automated algorithms to quan¬ 
tify the shape of a galaxy. One recent example of such a 
scheme can be found in Zehavi et al. (2002), who measured 
a “concentration parameter” for 30 000 galaxy images from 
the Sloan Digital Sky Survey, derived from the radii of differ¬ 
ent isophotes. Again, based upon cuts in the distribution of 
concentration parameter, early-types are found to be more 
clustered than late-types. 

In this paper, we employ a different method to classify 
galaxies, based upon a principal component analysis (PGA) 
of galaxy spectra, which is better suited to the 2dFGRS 
data (Madgwick et al. 2002). This technique has a number 
of attractive features. First, the PGA approach is completely 
objective and reproducible. An equivalent analysis can, for 
example, be applied readily to spectra produced by theo¬ 
retical models of galaxy formation or to spectra obtained 
in an independent redshift survey. Secondly, the PGA can 
be applied over the full magnitude range of the 2dFGRS, 
whenever the spectra are of sufficient signal to noise (see 
Section 2.2). For the 2dFGRS, the image quality is adequate 
to permit a visual determination of morphological type only 
for galaxies brighter than fej ~ 17, which comprise a mere 
5% of the spectroscopic sample. 

Two previous clustering studies have used spectral in¬ 
formation to select galaxy samples. Loveday, Tresse & Mad¬ 
dox (1999) grouped galaxies in the Stromlo-APM redshift 
survey into three classes based upon the equivalent width 
of either the Ha or OH lines, and found that galaxies 
with prominent emission lines display weaker clustering than 
more quiescent galaxies. Tegmark & Bromley (1999) mea¬ 
sured the relative bias between different spectral classes in 
the Las Campanas redshift survey (Shectman et al. 1996), 
using a classification derived from PCA analysis (Bromley et 
al. 1998), and also found that early spectral types are more 
strongly clustered than late spectral types. (See also Blan¬ 
ton 2000 for a revision of Tegmark & Bromley’s analysis, 
which takes into account the effect of errors in the survey 
selection function.) 

Here, we use the 2dFGRS survey to measure the depen¬ 
dence of galaxy clustering jointly on luminosity and spectral 


type, adding an extra dimension to the analysis carried out 
by Norberg et al. (2001). Previously, a pioneering study of 
bivariate galaxy clustering, in terms of luminosity and mor¬ 
phological type, was carried out using the Stromlo-APM 
redshift survey (Loveday et al. 1995). To place the analy¬ 
sis presented here in context, the samples that we consider 
cover a larg er volume and, despite being volume-limited (see 
Section 2.4), typically contain over an order of magnitude 


more galaxies than those available to Loveday et al. 

We give a brief overview of the 2dFGRS in Section ^ 
along with details of the spectral classification and an ex¬ 
planation of how the samples used in the clustering analy¬ 
sis were constructed. The estimation of the redshift space 
correlation function and its real space counterpartMhe pro¬ 


jected correlation function, are outlined in Section ^ A brief 
overview of the clustering of 2dFGRS galaxies in redshift 
space, selected by luminosity and spectral type, is given in 
Section a more detailed analysis of the redshift space clus¬ 
tering can be found in Hawkins et al. (2002). We present the 
main results of the paper in Section H and conclude in Sec¬ 
tion ^ 


2 THE DATA 

2.1 The 2dFGRS sample 

Detailed descriptions of the construction of the 2dFGRS and 
its properties are given by Colless et al. (2001). In summary, 
galaxies are selected down to a magnitude limit of bj « 
19.45 from the APM Galaxy Survey (Maddox et al. 1990a,b, 
1996, 2002). The sample considered in this paper consists of 
over 160 000 redshifts measured prior to May 2001. We focus 
our attention on the two large contiguous volumes of the 
survey, one centred on the Southern Galactic Pole (hereafter 
SGP) and the other close to the direction of the Northern 
Galactic Pole (NGP). 


2.2 Spectral classification of 2dFGRS galaxies 

The spectral properties of 2dFGRS galaxies are charac¬ 
terized using the principal component analysis (PGA) de¬ 
scribed by Madgwick et al. (2002). This analysis makes use 
of the spectral information in the rest-frame wavelength 
range 3700A to 6650A, thereby including all the major op¬ 
tical diagnostics between OH and Ha. For galaxies with 
2 > 0.15, sky absorption bands contaminate the Ha line. 
Since this can affect the stability of the classification, we re¬ 
strict our analysis to galaxies with 2 : < 0.15 following Madg¬ 
wick et al. ( 2002 ). 

The 2dFGRS spectra are classified by a single parame¬ 
ter, 77 , which is a linear combination of the hrst and second 
principal components. This combination has been chosen 
specihcally to isolate the relative strength of emission and 
absorption lines present in each galaxy’s spectrum, thereby 
providing a diagnostic which is robust to the instrumental 
uncertainties that affect the calibration of the continuum. 
Physically, this parameter is related to the star formation 
rate in a galaxy, as is apparent from the tight correlation 
of 77 with the equivalent width of Ha in emission line galax¬ 
ies (Bland-Hawthorn et al. 2002). In this paper, we divide 
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Figure 1. A comparison between the morphological classification 
of bright (bj < 17.0) APM galaxies by Loveday (1996) with the 
2dFGRS spectral classification, as quantified by the continuous 
variable rj (see text and Madgwick et al. 2002 for the definition). 
The morphological classification distinguishes between elliptical 
(Ell), lenticular (SO), spiral (Sp) and irregular (Irr) galaxies. All 
galaxies with both a morphological classification and a spectral 
classification are plotted. The Non-classified (Nc) class contains 
objects for which morphological classification was attempted but 
for which Loveday was unable to assign a type. The squares show 
the median value of r] for each morphological class defined by 
Loveday, and the error bars show the 10-90 percentile range. 



the 2dFGRS sample into two broad, distinct classes: galax¬ 
ies with spectra for which the PCA returns rj < —1.4, and 
which we refer to, for the sake of brevity, as early-type, and 
galaxies with rj > —1.4, which we call late-type. The distri¬ 
bution of v for 2dFGRS spectra displays a shoulder feature 
at this value (see Fig. 4 of Madgwick et al. 2002). 

The spectral type of a galaxy, as given by the value of 
rj, clearly depends upon its physical properties and is there¬ 
fore a useful and effective way in which to label galaxies. 
Nevertheless, it is still instructive to see how well, if at all, 
rj correlates with the morphogical type assigned in a sub¬ 
jective fashion from a galaxy image. Madgwick et al. (2002) 
show that there is a reasonable correspondence between rj 
and morphological type, using high signal-to-noise spectra 
and photometry taken from Kennicutt (1992); rj ~ —1.4 
approximately delineates the transition between early and 
late morphological types in bj. We revisit the comparison 
between classifications based on spectral and morphological 
types in Fig. this time using 2dFGRS spectra and UK 
Schmidt images. The horizontal axis shows the morpholog¬ 
ical type assigned to a subset of bright APM galaxies by 
Loveday (1996). Although there is a substantial amount of 
scatter in the rj values of spectra that lie within a given 
morphological class, it is reassuring to see that the median 
rj does correlate with morphological class. Moreover, the me¬ 
dian rj values match up well with the broad division that we 
employ to separate early and late types. Galaxies denoted 
“early-type” on the basis of their morphology have a median 
rj that is smaller than our fiducial value of »7 = —1.4 and vice- 
versa for late-types. In practice, for the samples analysed in 
this paper, the correspondence between morphological type 
and spectral class will be better than suggested by Fig. [U 
This is because the sample used in the comparison in Fig. ^ 
consists of nearby extended galaxies, and so the distribution 
of spectral types is distorted somewhat by aperture effects 
(see e.g. Kochanek, Pahre & Falco 2002; Madgwick et al. 
2002). This effect arises because the fibres used to collect 
the galaxy spectra are of finite size (subtending 2” on the 
sky). For this reason, when we measure the spectrum of a 
nearby galaxy it is possible that a disproportionate amount 
of light will be sampled from the bulge, thereby making the 
galaxy appear systematically redder or “earlier” in type. We 
find that this effect is only significant for the most nearby 
galaxies {z < 0.05) and should be completely negligible be¬ 
yond 2 ~ 0.1 (Madgwick et al. 2002). 


Figure 2. The histograms, plotted with Poisson error bars, show 
the success rate for assigning a spectral type to a targeted galaxy 
as a function of apparent magnitude. Two field completeness (cp, 
defined in text) ranges are shown, as indicated by the values at 
the top of each panel. The redshift completeness, i.e. the fraction 
of targeted galaxies for which a redshift is measured, is shown by 
the dotted lines. The spectral classification completeness, i.e. the 
fraction of galaxies with measured redshifts below z = 0.15 that 
have spectra of sufficient signal-to-noise to be used in the PCA, 
is shown by the dashed lines. (The dotted and dashed curves 
show parametric fits to the inferred redshift completeness and 
spectral classification completeness respectively.) The model for 
the spectral classification success rate, shown by the solid lines, 
is the product of the dotted and dashed lines in each panel, and 
is a good fit to the histogram in each case. 


2.3 Sample Selection 

In order to construct an optimal sample for the measure¬ 
ment of the two point correlation function, we select regions 
with high completeness in terms of measured redshifts, us¬ 
ing a redshift completeness mask for the 2dFGRS, similar to 
the one described in Colless et al. (2001; see also Norberg et 
al. 2002). Such a mask is required because of the tiling strat¬ 
egy adopted to make the best use of the allocated telescope 
time, along with the fact that the survey is not yet finished. 
An additional consideration is the success rate with which 
spectral types have been assigned to galaxies, which depends 
upon the signal-to-noise ratio of the galaxy spectrum. 

In Fig. 1^ we show histograms of the spectral classi¬ 
fication success rate for two different ranges of field com¬ 
pleteness, cf, which is defined as the ratio of the number 
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of measured redshifts in a given 2dF field to the number 
of targets. The spectral classification success rate has two 
contributions. The first of these is the redshift complete¬ 
ness, shown by the dotted curve. This incompleteness arises 
because we do not always succeed in measuring a redshift 
for a targeted object. The redshift incompleteness is nec¬ 
essarily small for the high completeness fields contributing 
to the histograms. The second contribution is the spectral 
classification completeness. Galaxies do not recieve a spec¬ 
tral classification when a redshift is measured with 2 < 0.15 
(and is therefore within the redshift range over which the 
PCA can be carried out), but the spectrum has too small 
a signal-to-noise ratio for the PCA to be applied success¬ 
fully (typically S/N < 10). The spectral classification suc¬ 
cess rate is given by the product of these two contributions. 
Our model for this effect, plotted as the solid curves in each 
panel of Fig. ^ is in good agreement with the success rate 
realised in the 2dFGRS, shown by the histograms. 

Rather than weight the data to compensate for a spec¬ 
tral classification success rate below 100%, we instead mod¬ 
ulate the number of unclustered or random points laid down 
in each field in the clustering analysis to take into account 
the varying success rate. We have conducted a number of 
tests in which we varied the completeness thresholds used, 
adopted different weighting schemes using samples of higher 
completeness, and we have also compared our results with 
those from Norberg et al. (2001), whose samples are not sub¬ 
ject to spectral classification incompleteness. The results of 
these tests confirm that our clustering measurements are 
robust to changes to the details of how we model the incom¬ 
pleteness; this is largely due to our practice of restricting the 
analysis to high completeness fields. Excluding areas below 
our relatively high sector completeness threshold (see Colless 
et al. 2001 for a definition), we estimate that the effective 
solid angle used in the SGP region is ^ 380 and in the 
NGP 250 0°. 

2.4 Constructing a volume-limited sample 

We analyse a series of volume-limited samples drawn from 
the 2dFGRS, following the strategy adopted by Norberg et 
al. (2001). The chief advantage of this approach is simplic¬ 
ity; the radial distribution of galaxies is uniform apart from 
modulations in space density due to clustering. Therefore, 
the complication of modelling the radial selection function in 
a flux-limited survey is avoided. This is particularly appeal¬ 
ing for the current analysis, as separate selection functions 
would be required for each class of spectral type studied, 
since Madgwick et al. (2002) have demonstrated that galax¬ 
ies with different spectral types have different luminosity 
functions. 

The disadvantage of using volume-limited samples is 
that a large fraction of galaxies in the flux-limited cata¬ 
logue do not satisify the selection criteria. As Norberg et 
al. (2001) point out, a volume-limited sample specified by 
a range in absolute magnitude has both a lower (zmin) and 
an upper redshift cut (zmax), because the flux-limited cata¬ 
logue has, in practice, bright and faint apparent magnitude 
limits. This seemingly profligate use of galaxy redshifts was 
a serious problem for previous generations of redshifts sur¬ 
veys. This is not, however, the case for the 2dFGRS, which 
contains sufficient galaxies to permit the construction of 


large volume-limited samples defined both by luminosity and 
spectral type. As we demonstrate in section ^ the volume- 
limited samples we analyse are large enough, both in terms 
of volume and number of galaxies, to give extremely robust 
clustering measurements. 

To construct a volume-limited sample, it is necessary 
to estimate the absolute magnitude that each galaxy would 
have at 2 = 0. This requires assumptions about the variation 
of galaxy luminosity with wavelength and with redshift, or 
equivalently, with cosmic time. We make use of the class 
dependent fc-corrections derived by Madgwick et al. (2002). 
The mean weighted fc-corrections are given by the following 
expressions: 


k{z) 

= 2.62 + 4.32^ (early — types) 

(1) 

k{z) 

= 1.5z + 2.1z^ (late — types) 

(2) 

k{z) 

= 1.9z + 2.7z^ (full sample). 

(3) 


These fc-corrections have the appeal that they are extracted 
directly in a self-consistent way from 2dFGRS spectra. How¬ 
ever, no account is taken of evolution in the galaxy spectrum. 
The explicit inclusion of evolution could lead to the am¬ 
biguous situation whereby a galaxy’s spectral type changes 
with redshift. We have checked that our results are, in fact, 
insensitive to the precise choice of fe-correction, comparing 
clustering results obtained with the spectral type dependent 
fc-corrections given above with those obtained when a global 
k + e-correction (i.e. making an explicit attempt to account 
for galaxy evolution, albeit in an average sense) is applied 
(as in Norberg et al. 2001). 

Since the fc-corrections are class dependent, the Zmin 
and 2max values corresponding to a given absolute magni¬ 
tude range are also slightly class dependent. Hence, the vol¬ 
umes defining the samples for two different spectral classes 
for the same bin in absolute magnitude will not coincide ex¬ 
actly. In addition to this subtle class dependent definition of 
the volumes, the values of 2min and 2max vary slightly with 
position on the sky. This is due to revisions made to the map 
of galactic extinction (Schlegel, Finkbeiner & Davis 1998) 
and to the CCD recalibration of APM plate zero-points since 
the definition of the original input catalogue. 

Finally, throughout the paper, we adopt an flo = 0.3, 
Ao = 0.7 cosmology to convert redshift into comoving dis¬ 
tance. The relative clustering strength of our samples is in¬ 
sensitive to this choice. 


3 ESTIMATING THE TWO-POINT 
CORRELATION FUNCTION 

The galaxy correlation function is estimated on a two dimen¬ 
sional grid of pair separation parallel (tt) and perpendicular 
(a) to the line-of-sight. To estimate the mean density of 
galaxy pairs, a catalogue of randomly positioned points is 
generated with the same angular distribution and the same 
values of 2min and 2max as the data, taking into account the 
completeness of the survey as a function of position on the 
sky, as described in Section 2.3, The correlation function is 
estimated using 


. DD RR 


( 4 ) 
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Figure 3. The spatial distribution of 2dFGRS galaxies in the SGP. The panels show the redshift and right ascension of galaxies in a 
three degree thick strip in declination for different magnitude ranges. To expand the scale of the plot, the redshift range shown has been 
restricted; note that the redshift scales are different in the two panels. Blue stars mark the locations of late-type galaxies and red circles 
the positions of early-type galaxies. The boxes show the two structures which are referred to in the text. 
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where DD, DR and RR are the numbers of data-data, data- 
random and random-random pairs respectively in each bin 
(Hamilton 1993). This estimator does not require an explicit 
estimate of the mean galaxy density. We have also cross¬ 
checked our results using the estimator proposed by Tandy 
& Szalay (1993): 


^LS 


DD - 2DR -b RR 


(5) 


where, this time DD, DR and RR are the suitably nor¬ 
malised numbers of data-data, data-random and random- 
random pairs. We find that the two estimators give the same 
results over the range of pair separations in which we are in¬ 
terested. 

The clustering pattern of galaxies is distorted when ra¬ 
dial positions are inferred from redshifts, as expected in the 
gravitational instability theory of structure formation (e.g. 
Kaiser 1987; Cole, Fisher & Weinberg 1994). Clear evidence 
for this effect is seen in the shape of the two point cor¬ 
relation function when plotted as ^(cr, tt), as demonstrated 
clearly for galaxies in the 2dFGRS by Peacock et al. (2001) 
and for groups of galaxies in the Zwicky catalogue by Padilla 
et al. (2001). After giving a brief flavour of the clustering of 
2dFGRS galaxies in redshift space in Section we focus 
our attention on clustering in real space in the remainder of 
the paper. The clustering signal in real space is inferred by 
integrating ^(cr, tt) in the tt direction (i.e. along the line of 
sight): 


H(a) 


a 


1 

- / C(cr, 7r)d7r. 
J —OO 


( 6 ) 


For the samples that we consider, the integral converges by 
pair separations of tt > 50/i“^Mpc. The projected corre¬ 
lation function can then be written as an integral over the 
spherically averaged real space correlation function, ^(r). 


S(cr) 


a 


2 

a 



rdr 

-a2)i/2’ 


(7) 


(Davis & Peebles 1983). If we assume that the real space 
correlation function is a power law (which is a fair ap¬ 
proximation for APM galaxies out to separations around 
r ~ 10h“^Mpc, see e.g. Baugh 1996), then Eq. can be 
written as 


E{a) 


(?)■ 


r(i/2)r([7-i]/2) 

r(7/2) 


( 0 )A,„, 


( 8 ) 


where r(a:) is the usual Gamma function, and we have used 
^(r) = (ro/r)^, where ro is the real space correlation length 
and 7 is equal to the slope of the projecte d cor relation func¬ 
tion H(ct) ja. As we demonstrate in Section 4.2, the projected 
correlation function is well described by a power law. 

We study a range of samples containing different num¬ 
bers of galaxies and covering different volumes of the Uni¬ 
verse. It is imperative to include sampling fluctuations when 
estimating the errors on the measured correlation function, 
to allow a meaningful comparison of the results obtained 
from different samples. This contribution to the errors has 
often been neglected in previous work. Following Norberg 
et al. (2001), we employ a sample of 22 mock 2dFGRS cat¬ 
alogues drawn from the AGDM Hubble Volume simulation 
(Evrard et al. 2002) to estimate the error bars on the mea¬ 
sured correlation functions. The construction of these mock 


catalogues is explained in Baugh et al. (2002, in preparation; 
see also Cole et al. 1998 and Norberg et al. 2002). These cat¬ 
alogues have the same selection criteria and the same clus¬ 
tering amplitude as measured for galaxies in the flux-limited 
2dFGRS. We have experimented with ensembles of mock 
catalogues constructed with different clustering strengths to 
ascertain how best to assign error bars when the measured 
clustering has a different amplitude from that of our fiducial 
sample of 22 2dFGRS mocks. We found that the error bars 
obtained directly by averaging over a test ensemble of mocks 
are reproduced most closely by using the scaled fractional 
rms scatter derived from the fiducial ensemble of 22 mocks, 
rather than by taking the absolute error. 


4 CLUSTERING IN REDSHIFT SPACE 


In this section we give a brief overview of the clustering of 
2dFGRS galaxies in redshift space, for sampl es s elected by 
luminosity and spectral type. First, in Section 4.1, we give a 


qualitative impression of the clustering differences by plot¬ 
ting the spatial distribution of galaxies in volume-limited 
samples. Then we quantify these differences by measuring 
the spherically averaged correlation function, ^(s). A more 
comprehensive analysis of the clustering of 2dFGRS galaxies 
in redshift space will be presented by Hawkins et al. (2002). 


4.1 Spatial distribution of 2dFGRS galaxies 

It is instructive to gain a visual impression of the spa¬ 
tial distribution of 2dFGRS galaxies before interpreting the 
measured correlation functions. In Fig. ^ we show the spa¬ 
tial distribution of galaxies in two ranges of absolute mag¬ 
nitude: in panel (a) we show a sample of faint galaxies 
( — 18.0 > Mbj—SlogjQ h > —19.0) and in panel (b) a sam¬ 
ple of bright galaxies (—20.0 > Mhj — Slog^Q h > —21.0). 
Within each panel, early and late type galaxies, as distin¬ 
guished by their spectral types, are plotted with different 
symbols; the positions of early-types are indicated by circles 
and the late-types are marked by stars. For clarity, we show 
only a three degree slice in declination cut from the SGP re¬ 
gion and we have sparsely sampled the galaxies, so that the 
space densities of the two spectral classes are equal. In order 
to expand the scale of the plot, the range of redshifts shown 
is restricted, taking a subset of the full volume-limited sam¬ 
ple in each case. (Note also that the redshift ranges differ 
between the two panels.) 

A hierarchy of structures is readily apparent in these 
plots, ranging from isolated objects, to groups of a handful 
of galaxies and on through to rich clusters containing over 
a hundred members. It is interesting to see how structures 
are traced by galaxies in the different luminosity bins by 
comparing common structures between the two panels. For 
example, the prominent structure (possibly a supercluster 
of galaxies) seen at a ~ 0^ and a ~ 0.061 is clearly visible 
in both panels. The same is true for the overdensity seen at 
a ~ 03‘'15' at z ~ 0.068. 

This is the first time that a large enough survey has 
been available, both in terms of the volume spanned and 
the number of measured redshifts, to allow a comparison 
of the clustering of galaxies of different spectral types and 
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luminosities in representative volume-limited samples, with¬ 
out the complication of the strong radial gradient in number 
density seen in flux-limited samples. 

It is apparent from a comparison of the distribution 
of the different spectral types in Fig. |^(a), that the faint 
early-type galaxies tend to be grouped into structures on 
small scales whereas the faint late-types are more spread out. 
One would therefore anticipate that the early-types should 
have a stronger clustering amplitude tha n th e late-types, an 
expectation that is borne out in Section 4.2. 

In Fig. f 


the distinction between the distribution of 
the spectral types is less apparent. This is partly due to the 
greater importance of projection effects in the declination 
direction, as the cone extends to a greater redshift than in 
Fig. ^(a). However, close examination of the largest struc¬ 
tures suggests that early-types are more abundant in them 
than late-types, again implying a stronger clustering ampli¬ 
tude. 


4.2 ^(s) as function of luminosity and spectral 

type 

In Fig. ^ we show the spherically averaged redshift space 
correlation function, ^(s), as a function of luminosity and 
spectral type. Results are shown for samples selected in bins 
of width one magnitude, as indicated by the legend in each 
panel. The top panel shows the correlation functions of all 
galaxies that have been assigned a spectral type, the mid¬ 
dle panel shows results for galaxies classified as early-types 
(rj < —1.4) and the bottom panel shows results for late-types 
(rj > —1.4). Note that, at present, there are insufficient num¬ 
bers of late-type galaxies to permit a reliable measurement 
of the correlation function for the brightest magnitude bin 
-21.0 > Mtj - Slogio h > -22.0. 

Several deductions can be made immediately from 
Fig. ^ In all cases, the redshift space correlation func¬ 
tion is well described by a power-law only over a fairly 
limited range of scales. The correlation functions of early- 
type galaxies are somewhat steeper than those of late- 
types. However, the main difference is that the early-type 
galaxies have a stronger clustering amplitude than the late- 
type galaxies. The correlation length, defined here as the 
pair separation for which ^(so) = 1, varies for early- 

types from So = 7.1 ± 0.7/i“^Mpc for galaxies with ab¬ 
solute magnitudes around M(,j — SlogjQ h ^ —19.5 to 
So = 8.9 ± 0.7/i~^Mpc for the brightest sample with 
—21.0 > Mbj—Slogj^g h > —22.0. The faintest early-types, 
with magnitudes—17.5 > M^j-51og^Q h > —18.5, display 
a clustering amplitude that is similar to that of the bright¬ 
est early-types. However, the measurement of the correlation 
function for this faint sample is relatively noisy, as the vol¬ 
ume in which galaxies are selected is small compared with 
the volumes used for brighter samples. The late-type galax¬ 
ies show, by contrast, little change in clustering amplitude 
with increasing luminosity, with a redshift space correlation 
length of So = 5.6 ± 0.6/i“^Mpc. Only a slight steepen¬ 
ing of the redshift space correlation function is apparent 
with increasing luminosity, until the brightest sample, which 
displays a modest increase in the redshift space correlation 
length. The correlation lengths of all our samples of early- 
type galaxies are larger than those of late-type galaxies. 
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Figure 4. The spherically averaged redshift space correlation 
function of galaxies in disjoint absolute magnitude bins, as in¬ 
dicated by the key in each panel. The panels show the results 
for different samples: the top panel shows the correlation func¬ 
tions for all galaxies that have been assigned a spectral type, 
the middle panel show the clustering of galaxies with 77 < —1.4 
and the lower panel shows ^(s) for galaxies with rj > —1.4. The 
error bars are obtained using 2dFGRS mock catalogues, as de¬ 
scribed in the text. For clarity, error bars are plotted only on the 
— 18.5 > Affcj — 51og;LQ h > —19.5 sample curve and for the 
brightest sample in each panel. 


5 CLUSTERING IN REAL SPACE 

5.1 Robustness of clustering results 

The approach adopted to study the dependence of galaxy 
clustering on luminosity relies upon being able to compare 
correlation functions measured in different volumes. It is im¬ 
portant to ensure that there are no systematic effects, such 
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Table 1. Properties of the combined NGP SGP volume-limited samples for all galaxies with a spectral type. Column 1 lists the 
absolute magnitude range which defines the volume-limited samples. Columns 2 and 3 list the median magnitude and the number of 
galaxies in the sample. Columns 4 and 5 give the redshift limits of the sample. Columns 6 and 7 list the best fitting correlation length, 
ro, and power-law slope, 7 , of the real space correlation function, and column 8 gives ^( 7 ), as defined by Eq. 8 . 


Mag. range 

Mjj -51ogiofi 

Median magnitude 
Mjj - 51ogiofi 

-Ngal 

■^min 

• 2 max 

ro 

(/i-^Mpc) 

7 

±( 7 ) 

-17.5 

- 18.5 

-17.98 

8510 

0.0164 

0.0724 

5.19 ±0.95 

1.68 ± 0.12 

4.14 

-18.0 

- 19.0 

-18.46 

13795 

0.0204 

0.0886 

4.36 ±0.89 

1.83 ± 0.10 

3.58 

-18.5 

- 19.5 

-18.93 

19207 

0.0255 

0.1077 

4.65 ±0.61 

1.80 ± 0.08 

3.68 

-19.0 

- 20.0 

-19.40 

24675 

0.0317 

0.1302 

4.93 ±0.50 

1.79 ± 0.10 

3.71 

-19.5 

- 20.5 

-19.85 

22555 

0.0394 

0.1500 

4.89 ±0.48 

1.79 ± 0.05 

3.71 

- 20.0 

- 21.0 

-20.30 

10399 

0.0487 

0.1500 

5.37 ±0.61 

1.78 ± 0.11 

3.75 

-20.5 

- 21.5 

-20.74 

3423 

0.0602 

0.1500 

6.57 ±0.83 

1.83 ± 0.23 

3.58 

- 21.0 

- 22.0 

-21.19 

751 

0.0739 

0.1500 

8.47 ±0.97 

1.80 ± 0.29 

3.68 


Table 2. Properties of the combined NGP & SGP volume-limited samples of early-type galaxies. See Table 1 for column definitions. 


Mag. range Median magnitude -^gal ^max 7 ^( 7 ) 

Mij-Slogigh M6j-51ogiQh (ft-iMpc) 


-17.5 

- 18.5 

-18.05 

1909 

0.0163 

0.0707 

8.33 ± 1.82 

1.87 ±0.23 

3.46 

-18.0 

- 19.0 

-18.53 

3717 

0.020.3 

0.0861 

6.28 ± 1.46 

1.98 ± 0.11 

3.19 

-18.5 

- 19.5 

-18.98 

6405 

0.0253 

0.1041 

5.92 ± 1.00 

1.83 ± 0.10 

3.58 

-19.0 

- 20.0 

-19.44 

10135 

0.0314 

0.1249 

5.71 ±0.57 

1.87 ±0.09 

3.46 

-19.5 

- 20.5 

-19.89 

11346 

0.0.388 

0.1486 

5.66 ±0.56 

1.87 ±0.09 

3.46 

- 20.0 

- 21.0 

-20.33 

6434 

0.0480 

0.1500 

6.10 ±0.72 

1.80 ± 0.12 

3.68 

-20.5 

- 21.5 

-20.77 

2587 

0.0590 

0.1500 

7.60 ± 1.02 

1.87 ±0.26 

3.46 

- 21.0 

- 22.0 

- 21.21 

686 

0.0722 

0.1500 

9.74 ± 1.16 

1.95 ± 0.37 

3.26 


Table 3. Properties of the combined NGP & SGP volume-limited samples of late-type galaxies. See Table 1 for column definitions. Note 
that the brightest sample listed contains too few galaxies to permit a reliable measurement of the projected correlation function. 


Mag. range 

Mbj -51ogioh 

Median magnitude 

Mjj - 51ogiQ h 

■Ngal 

•^min 

•^max 

ro 

(/i-^Mpc) 

7 

A(7) 

-17.5 

- 18.5 

-17.96 

6674 

0.0164 

0.0734 

4.27 ±0.81 

1.65 ± 0.12 

4.29 

-18.0 

- 19.0 

-18.44 

9992 

0.0205 

0.0901 

3.71 ±0.77 

1.76 ± 0.11 

3.82 

-18.5 

- 19.5 

-18.90 

12619 

0.0256 

0.1099 

4.17 ±0.64 

1.79 ± 0.10 

3.71 

-19.0 

- 20.0 

-19.37 

14420 

0.0319 

0.1333 

4.45 ±0.47 

1.76 ± 0.09 

3.82 

-19.5 

- 20.5 

-19.82 

11122 

0.0397 

0.1500 

4.59 ±0.44 

1.76 ± 0.07 

3.82 

- 20.0 

- 21.0 

-20.26 

4300 

0.0492 

0.1500 

5.52 ±0.88 

1.87 ±0.1.3 

3.46 

-20.5 

- 21.5 

-20.71 

1118 

0.0608 

0.1500 

6.33 ± 1.01 

2.01 ± 0.29 

3.12 

- 21.0 

- 22.0 

-21.17 

198 

0.0749 

0.1500 

— 

- 

— 


as significant sampling fluctuations, that could undermine 
such an analysis. In Norberg et al. (2001) we demonstrated 
the robustness of this approach in two ways. First, we con¬ 
structed a volume-limited sample defined using a broad mag¬ 
nitude range, that could be divided into co-spatial subsam¬ 
ples of galaxies in different luminosity bins, i.e. subsamples 
within the same volume and therefore subject to the same 
large-scale structure fluctuations. A clear increase in clus¬ 
tering amplitude was found for the brightest galaxies in the 
volume, establishing the dependence of clustering on galaxy 
luminosity (see Fig. la of Norberg et al. 2001). Secondly, 
we demonstrated that measuring the correlation function of 
galaxies in a fixed luminosity bin, but using samples taken 
from different volumes, gave consistent results (see Fig. lb 
of Norberg et al. 2001). 

In this section we repeat these tests. The motivation 
for this exercise is that the samples considered here contain 
fewer galaxies than those used by Norberg et al. (2001) as 


only galaxies with z < 0.15 are suitable for PCA spectral 
typing, and because the samples are more dilute as they 
have been selected on the basis of spectral type as well as 
luminosity. In Fig. |^a) we plot the projected correlation 
function of late-type galaxies in a fixed absolute magnitude 
bin ( — 19.0 > M(,j — Slogj^Q h > —20.0), but measured for 
samples taken from volumes defined by different ^min and 
2max. The clustering results are in excellent agreement with 
one another. In Fig. |^(b), we compare the projected cor¬ 
relation function of late-type galaxies in different absolute 
magnitude ranges but occupying the same volume. A clear 
difference in the clustering amplitude is seen. We have also 
performed these tests for early-type galaxies and arrived at 
similar conclusions. 

As an additional test, we also show in Fig. ^ the cor¬ 
relation function measured in what we refer to as the opti¬ 
mal sample for a given magnitude bin. The optimal sample 
contains the maximum number of galaxies for the specified 
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Figure 5. (a) The projected correlation function of late-type 

galaxies in a fixed absolute magnitude bin taken from different, 
almost independent volumes. We show the correlation function of 
galaxies with —19.0 > Mjj — Slogj^Q h > —20.0 taken from 
volumes defined by —18.0 > Mi,^ — SlogjQ h > —20.0 and 

— 19.0 > Mjj-51ogjg h > —21.0 (both shown by heavy dashed 
lines). The thin solid line shows the estimate from the optimal 
sample for the —19.0 > Mjj — hlogj^g h > —20.0 magnitude 
bin. The different measurements are in almost perfect agreement, 
(b) The projected correlation function measured for late-type 
galaxies in two different absolute magnitude bins taken from the 
same volume. The volume is defined by the magnitude range 

— 19.0 > Mjj — Slogj^Q h > —21.0. Within a fixed volume, there 
is clear evidence for an increase (albeit small) in the clustering 
amplitude with luminosity. The two thin solid lines show esti¬ 
mates obtained from the corresponding optimal samples for the 
stated magnitude bins. In both panels the error bars come from 
the analysis of mock 2dFGRS catalogues. 


magnitude bin. The correlation functions of galaxies in op¬ 
timal samples are shown by thin solid lines in both panels 
and are in excellent agreement with the other measurements 
shown. 


5.2 Projected correlation function 

Fig. ^ shows how the real space clustering of galaxies of dif¬ 
ferent spectral type depends on luminosity. We use the opti¬ 
mal sample for each magnitude bin, i.e. the volume-limited 
sample with the maximum possible number of galaxies, the 
properties of which are listed in Tables (early & late types 
together), ^ (early-types only) and ^ (late-types only). 


The top panel of Fig. ^ confirms the results found by 
Norberg et al. (2001), namely that the clustering strength of 
the full sample increases slowly with increasing luminosity 
for galaxies fainter than M* , and then shows a clear, strong 
increase for galaxies brighter than M* . (We take M* to be 
Mt,j — Slogj^Q h ~ —19.7, following Folkes et al. 1999 and 
Norberg et al. 2002.) Furthermore, the projected correla¬ 
tion functions are well described by a power law with a slope 
that is independent of luminosity. The middle panel of Fig. ^ 
shows the projected correlation function of early-type galax¬ 
ies for different absolute magnitude ranges. The clustering 
amplitude displays a non-monotonic behaviour, with the 
faintest sample having almost the same clustering strength 
as the brightest sample. The significance of this result for 
the faintest galaxies will be discussed further in the next 
section. Early-type galaxies with Mbj — Slogj^g ^ — “19.5, 
display weaker clustering than the faint and bright samples. 
The bottom panel of Fig. ^ shows the real space clustering 
of late-type galaxies as a function of luminosity. In this case, 
the trend of clustering strength with luminosity is much 
simpler. There is an increase in clustering amplitude with 
luminosity, and also some evidence that the projected cor¬ 
relation function of the brightest subset is steeper than that 
of the other late-type samples. In general, for the luminosity 
ranges for which a comparison can be made, the clustering 
strength of early-type galaxies is always stronger than that 
of late-types. 

The comparison of the correlation functions of differ¬ 
ent samples is made simpler if we divide the curves plot¬ 
ted in Fig. ^ by a fiducial correlation function. As a ref¬ 
erence sample we choose all galaxies that have been as¬ 
signed a spectral type, with absolute magnitudes in the 
range —19.5 > M(,j — 51ogj^n h > —20.5 (the short-dashed 
line in the top panel of Fig. In Fig. we plot the ratio 
of the correlation functions shown in the panels of Fig. ^ to 
the reference correlation function, with error bars obtained 
from the mock catalogues. The trends reported above for the 
variation of clustering strength with luminosity and spec¬ 
tral type are now clearly visible (see Fig. |^, particularly 
the difference in clustering amplitude between early-types 
and late-types. In the upper and lower panels, the ratios of 
correlation functions are essentially independent of a, indi¬ 
cating that a single power-law slope is a good description 
over the range of scales plotted. The one exception is the 
brightest sample of late-type galaxies, which shows some 
evidence for a steeper power-law. In the middle panel, the 
ratios for early-type galaxies show tentative evidence for a 
slight steepening of the correlation function at small pair 
separations, a < 2h~^ Mpc, which is most pronounced for 
the brightest sample. 


5.3 Real space correlation length 

In the previous subsection, we demonstrated that the pro¬ 
jected correlation functions of galaxies in the 2dFGRS have 
a power-law form with a slope that varies little as the sample 
selection is changed, particularly for pair separations in the 
range 2.0 < a/{h~^ Mpc) < 15.0. To summarize the trends 
in clustering strength found when varying the spectral type 
and luminosity of the sample, we fit a power-law over this 
range of scales. We follow the approach, based on Eq. used 
by Norberg et al. (2001), who performed a minimisation 
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Figure 6. The projected galaxy correlation function for sam¬ 
ples of different spectral type, split into one absolute magnitude 
wide bins. The top panel shows the correlation function measured 
for all galaxies with an assigned spectral type. The middle panel 
shows correlation functions for early-types and the bottom panel 
shows the results for late-types. The absolute magnitude range 
of each sample is indicated in the legend on each panel. The 
error bars are derived from the 2dFGRS mock catalogues and, 
for clarity, are only plotted on the correlation functions of the 
— 18.5 > — 5 logj^Q h > —19.5 sample and of the brightest 

samples in each panel. 


to extract the best fitting values of the parameters in the 
power-law model for the real space correlation function: the 
correlation length, ro, and the power law slope, 7 . As pointed 
out by Norberg et al. (2001), a simple approach does not 
give reliable estimates of the errors on the fitted parameters 
because of the correlation between the estimates of the cor¬ 
relation function at differing pair separations. We therefore 



-0.5 0 0.5 1 

logio(cr/h"^Mpc) 


Figure 7. The ratio of the projected correlation function mea¬ 
sured for galaxies selected by luminosity and spectral type to the 
projected correlation function of a reference sample. The refer¬ 
ence sample consists of all galaxies with an assigned spectral type 
that lie within the magnitude range —19.5 > Mf,^ — hlog^Q h > 
—20.5. The top panel shows the ratios for different luminosity 
bins for all galaxies with a spectral type, the middle panel shows 
the ratios obtained for early-types and the bottom panel shows 
the ratios for late-types. The same line styles plotted in Fig. ^ 
are used to indicate different luminosities. The error bars, which 
take into account the correlation between the samples (but not 
between the bins), are from the mock 2dFGRS catalogues, and 
for clarity, are only plotted on two curves in each panel: the ratio 
for the sample with —18.5 > M(,j — 51og]^g h > —19.5 and the 
ratio for the brightest sample in each panel. 


use the mock 2dFGRS catalogues to estimate the errors on 
the fitted parameters in the following manner. The best fit¬ 
ting values of rg and 7 are found for each mock individually, 
using the analysis. The estimated error is then taken to 
be the fractional rms scatter in the fitted parameters over 
the ensemble of mock catalogues. The best fitting parame¬ 
ters for each 2dFGRS sample are listed in Tables llJ (early & 
late types together), (early-types only) and (late-types 
only). The results for the correlation length are plotted in 
the top panel of Fig. 

The correlation lengths estimated for the full sample 
with assigned spectral types (shown by the open squares in 
Fig. 1^ are in excellent agreement with the results of Nor¬ 
berg et al. (2001) (shown by the filled circles). The bright 
samples constructed by Norberg et al. (2001) have amax val¬ 
ues in excess of the limit of Zmax= 0.15 enforced upon the 
samples analysed in this paper by the PCA. The bright sam¬ 
ples used in this paper therefore cover smaller volumes than 
those used by Norberg et al. and so the error bars are sub- 
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Figure 8. Top panel: the real space correlation length, tq, as 
a function of the absolute magnitude of the sample, for galaxies 
of different spectral type. The stars show the tq values fitted to 
the projected correlation function of late-type galaxies and the 
open circles show the ro values for early-types. The squares show 
ro for the full sample with spectral types. The latter results are 
in excellent agreement with those obtained from the larger sam¬ 
ple analysed by Norberg et al. (2001), which are plotted as filled 
circles. The horizontal bars on the filled circles show the mag¬ 
nitude range used to define the volume-limited samples. Lower 
panel: the difference Ar^ = rQ(Mbj) — for each spec¬ 

tral type. This quantity gives the significance of the variation of 
the correlation length with luminosity for each spectral type with 
respect to a reference sample, chosen to be the sample defined 
by —19.5 > Mjj — hlogj^Q h > —20.5. The error bars take into 
account the correlation between the two samples. 


stantially larger. A cursory inspection of Fig. ^ would give 
the misleading impression that we find weaker evidence for 
an increase in correlation length with luminosity. It is im¬ 
portant to examine this plot in conjunction with Tables ^ 
to 1^ which reveals that there is significant overlap in the 
volumes defined by the four brightest magnitude slices, be¬ 
cause of the common Zmax= 0.15 limit. 

In this case, the error bars inferred directly from the 
mocks do not take into account that our samples are corre¬ 


lated. The errors fully incorporate cosmic variance, i.e. the 
variance in clustering signal expected when sampling a given 
volume placed at different, independent locations in the Uni¬ 
verse. The volumes containing the four brightest samples 
listed in Table contain long-wavelength fluctuations in 
common and so clustering measurements from these differ¬ 
ent volumes are subject to a certain degree of coherency. 
The clearest way to show this is by calculating the differ¬ 
ence between the correlation lengths fitted to two samples. 
The error on the difference, derived using the mock cata¬ 
logues, has two components: the first comes from adding the 
individual errors in quadrature; the second is from the cor¬ 
relation of the samples. For correlated samples, this second 
term is negative and, therefore lowers the estimated error 
on the difference. This is precisely what is seen in the lower 
panel of Fig. ^ where we plot Ar^ — rQ(Mbj) — rQ(M[j®j) 
for each spectral type as a function of absolute magnitude, 
with error bars taking into account the correlation of the 
samples. For all samples brighter than M* , the increase in 
clustering length with luminosity is clear. 

There is a suggestion, in the top panel of Fig. of a non¬ 
monotonic dependence of the correlation length on luminos¬ 
ity for early-type galaxies. The evidence for this behaviour is 
less apparent on the A panel, where a difference in rg for 
the faintest sample is seen at less than the 2cr level. These 
volumes are small compared to those defining the brighter 
samples. Furthermore, when analyzing the faintest sample 
for SGP and NGP separately, our results are somewhat sen¬ 
sitive to the presence of single structures in each region (at 
a ~ O*' and z ~ 0.061 for the SGP, as shown in Fig. ^(a)). 
Thus we conclude that the upturn at faint magnitudes in 
the correlation length of early-types is not significant. 

The projected correlation function of early-type galax¬ 
ies brighter than M* is well fitted by a power-law real 
space correlation function, with a virtually constant slope 
of 7 ~ 1.87 and a correlation length which increases with 
luminosity, from ro = 5.7 ± 0.6/i“^Mpc for M* galaxies 
to ro = 9.7 ± 1.2/i“^Mpc for brighter galaxies (Mtj — 
51og]^Q h ~ —21.2). This represents an increase in clustering 
strength by a factor of 2.7, as seen in Fig. |^. The projected 
correlation functions of late-type galaxies are also consistent 
with a power-law in real space, with an essentially constant 
slope. There is a very weak trend for 7 to increase with lu¬ 
minosity, although at little more than the Ict level. Ignoring 
this effect, the fitted slope of the late-type correlation func¬ 
tion is 7 ~ 1.76 . The correlation length increases with lumi¬ 
nosity from a value of ro = 3.7 ± 0.8/i“^Mpc for faint galax¬ 
ies {Mbj — blogj^Q h ~ —18.4) to ro = 6.3 ± 1.0/i“^Mpc 
for bright galaxies (M(,j — blog^o h ~ —20.7), a factor 
of 2.5 increase in clustering strength. It should be pos¬ 
sible to extend the analysis for late-type galaxies beyond 
Mb] — 51og]^o — —21 when the 2dFGRS is complete. 

The top panel of Fig. ^ confirms our earlier conclu¬ 
sion that the clustering of early-type galaxies is stronger 
than that of late-type galaxies. At M* , early-types typically 
have a real space clustering amplitude that is 1.5 — 1.7 times 
greater than that of late-types. 
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Figure 9. The fraction of galaxies in the two broad spectral 
classes, early-type and late-type, as a function of absolute magni¬ 
tude. The fractions are derived from the volume-limited samples 
listed in Tables 1, 2 and 3. The error bars show the Poisson errors 
on the fractions. 


6 DISCUSSION AND CONCLUSIONS 

We have used the 2dFGRS to study the dependence of clus¬ 
tering on spectral type for samples spanning a factor of 
twenty in galaxy luminosity. The only previous attempt at 
a bivariate luminosity-morphology/spectral type analysis of 
galaxy clustering was performed by Loveday et al. (1995) 
using the Stromlo-APM redshift survey. They were able to 
probe only a relatively narrow range in luminosity around 
L*, which is more readily apparent if one considers the me¬ 
dian magnitude of each of their magnitude bins (see Fig. 3b 
of Norberg et al. 2001). The scatter between spectral and 
morphological types illustrated in Fig. precludes a more 
detailed comparison of our results with those of earlier stud¬ 
ies, based on morphological classifications. 

In Norberg et al. (2001), we used the 2dFGRS to make a 
precise measurement of the dependence of galaxy clustering 
on luminosity. The clustering amplitude was found to scale 
linearly with luminosity. One of the aims of the present pa¬ 
per is to identify the phenomena that drive this relation. In 
particular, there are two distinct hypotheses that we wish to 
test. The first is that there is a general trend for clustering 
strength to increase with luminosity, regardless of the spec¬ 
tral type of the galaxy. The second is that different types of 
galaxies have different clustering strengths, which may vary 
relatively little with luminosity, but a variation of the mix 
of galaxy types with luminosity results in a dependence of 
the clustering strength on luminosity. 

Madgwick et al. (2002) estimated the luminosity func¬ 
tion of 2dFGRS galaxies for different spectral classes, and 
found that, in going from early-type to late-types, the slope 
of the faint end of the luminosity function becomes steeper 
while the characteristic luminosity becomes fainter. Another 
representation of the variation of the luminosity function 
with spectral class is shown in Fig. ^ where we plot the frac¬ 
tion of early and late type galaxies in absolute magnitude 
bins. The plotted fractions are derived from the volume- 


Figure 10. The relative bias (on the scale of r = 4.89/r“^Mpc) 
of the different spectral classes as a function of luminosity, as in¬ 
dicated by the key. The definition of relative bias is given in the 
text. The reference sample covers the absolute magnitude range 
— 19.5 > — 51og]^g h > —20.5. The fiducial luminosity, L*, 

is taken to be M(,j — Slogj^g h = —19.7. The solid line shows the 
fit to the results of Norberg et al. (2001), b/b* = 0.85 + 0.15L/L*, 
whereas the dashed line shows the fit to the data analyzed here. 
All the error bars plotted take into account the correlation be¬ 
tween the various samples. 


limited samples listed in Tables 1, 2 and 3. The mix of spec¬ 
tral types changes dramatically with luminosity; faint sam¬ 
ples are dominated by late-types, whereas early-types are 
the most common galaxies in bright samples. Similar trends 
were found for galaxies labelled by morphological type in 
the SSRS2 survey by Marzke et al. (1998). 

We find that the change in the mix of spectral types 
with luminosity is not the main cause for the increase in 
the clustering strength of the full sample with luminosity. 
To support this assertion, we plot in Fig. the variation 
of clustering strength with luminosity normalized, for each 
spectral class, to the clustering strength of a fiducial sample 
of M* galaxies, i. e. the sample which covers the magnitude 
range —19.5 > Mj,j — 51og^g h > —20.5. For a galaxy 
sample with best fitting correlation function parameters Vq 
and 7 *, we define the relative bias with respect to the M* 
sample of the same type by 



{rhT 




(9) 


where ro and 7 are the best fitting power-law parameters 
for the fiducial sample. In Fig. E we plot the relative bias 
evaluated at a fixed scale, r = 4.89h“^Mpc, which is the 
correlation length of the reference sample for all 77 -classified 
galaxies. A scale dependence in Eq. arises if the slopes 
of the real space correlation functions are different for the 
galaxy samples being compared. In practice, the term 
is close to unity for the samples considered. The dashed line 
shows a fit to the bias relation defined by the open symbols. 
The solid line shows the effective bias relation obtained by 
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Norberg et al. (2001), which is dehned in a slightly different 
way to the effective bias computed here. 

From Fig. ^ we see that the trend of increasing clus¬ 
tering strength with luminosity in both spectral classes is 
very similar for galaxies brighter than L > 0.5L*. At the 
brightest luminosity, corresponding to ~ 4L*, the cluster¬ 
ing amplitude is a factor of 2 — 2.5 times greater than at L*. 
This increase is much larger than the offset in the relative 
bias factors of early and late types at any given luminosity. 
We conclude that the change in correlation length with ab¬ 
solute magnitude found by Norberg et al. (2001) is primarily 
a luminosity effect rather than a reflection of the change in 
the mix of spectral types with luminosity. 

Benson et al. (2001) showed that a dependence of clus¬ 
tering strength on luminosity is expected in hierarchical 
clustering cold dark matter universes because of the pref¬ 
erential formation of the brightest galaxies in the most mas¬ 
sive, strongly clustered dark halos. The close connection be¬ 
tween the spectral characteristics of galaxies and their clus¬ 
tering properties discussed in this paper provides further 
evidence that the galaxy type is also related to the mass of 
the halo in which galaxy forms. 
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Table 1. Properties of the combined NGP & SGP volume-limited samples for all galaxies with a spectral type. Column 1 lists the 
absolute magnitude range which defines the volume-limited samples. Columns 2 and 3 list the median magnitude and the number of 
galaxies in the sample. Columns 4 and 5 give the redshift limits of the sample. Columns 6 and 7 list the best fitting correlation length, 
ro, and power-law slope, 7 , of the real space correlation function, and column 8 gives ^( 7 ), as defined by Eq. 0 


Mag. range 
Mhj - Slog^Q h 

Median magnitude 
Mhj - 5 logjo h 

^ga.1 

-^min 

-^max 

ro 

(/i“^Mpc) 

7 

Ml) 

-17.5 

- 18.5 

-17.98 

8510 

0.0164 

0.0724 

5.19 ± 0.95 

1.68 ±0.12 

4.14 

-18.0 

- 19.0 

-18.46 

13795 

0.0204 

0.0886 

4.36 ±0.89 

1.83 ±0.10 

3.58 

-18.5 

- 19.5 

-18.93 

19207 

0.0255 

0.1077 

4.65 ±0.61 

1.80 ±0.08 

3.68 

-19.0 

- 20.0 

-19.40 

24675 

0.0317 

0.1302 

4.93 ±0.50 

1.79 ±0.10 

3.71 

-19.5 

- 20.5 

-19.85 

22555 

0.0394 

0.1500 

4.89 ± 0.48 

1.79 ±0.05 

3.71 

-20.0 

- 21.0 

-20.30 

10399 

0.0487 

0.1500 

5.37 ±0.61 

1.78 ±0.11 

3.75 

-20.5 

- 21.5 

-20.74 

3423 

0.0602 

0.1500 

6.57 ±0.83 

1.83 ±0.23 

3.58 

-21.0 

- 22.0 

-21.19 

751 

0.0739 

0.1500 

8.47 ±0.97 

1.80 ±0.29 

3.68 


Table 2. Properties of the combined NGP Sz SGP volume-limited samples of early-type galaxies. See Tabled for column definitions. 


Mag. range Median magnitude Ngai -^min -^max ^0 7 ^( 7 ) 

Mt,j-5logiQh M6j-51ogig/i (fe'^Mpc) 


-17.5 

- 18.5 

-18.05 

1909 

0.0163 

0.0707 

8.33 ± 1.82 

1.87 ±0.23 

3.46 

-18.0 

- 19.0 

-18.53 

3717 

0.0203 

0.0861 

6.28 ± 1.46 

1.98 ±0.11 

3.19 

-18.5 

- 19.5 

-18.98 

6405 

0.0253 

0.1041 

5.92 ± 1.00 

1.83 ±0.10 

3.58 

-19.0 

- 20.0 

-19.44 

10135 

0.0314 

0.1249 

5.71 ±0.57 

1.87 ±0.09 

3.46 

-19.5 

- 20.5 

-19.89 

11346 

0.0388 

0.1486 

5.66 ±0.56 

1.87 ±0.09 

3.46 

- 20.0 

- 21.0 

-20.33 

6434 

0.0480 

0.1500 

6.10 ±0.72 

1.80 ± 0.12 

3.68 

-20.5 

- 21.5 

-20.77 

2587 

0.0590 

0.1500 

7.60 ± 1.02 

1.87 ±0.26 

3.46 

- 21.0 

- 22.0 

- 21.21 

686 

0.0722 

0.1500 

9.74 ± 1.16 

1.95 ±0.37 

3.26 


Table 3. Properties of the combined NGP Sz SGP volume-limited samples of late-type galaxies. See Table for column definitions. 
Note that the brightest sample listed contains too few galaxies to permit a reliable measurement of the projected correlation function. 


Mag. range 
Mj,j - 51ogio h 

Median magnitude 
Mhj - 5 logio h 

^ga.1 

^min 

-^max 

?'o 

(fe~^Mpc) 

7 

A ( 7 ) 

-17.5 

- 18.5 

-17.96 

6674 

0.0164 

0.0734 

4.27 ±0.81 

1.65 ±0.12 

4.29 

-18.0 

- 19.0 

-18.44 

9992 

0.0205 

0.0901 

3.71 ±0.77 

1.76 ±0.11 

3.82 

-18.5 

- 19.5 

-18.90 

12619 

0.0256 

0.1099 

4.17 ±0.64 

1.79 ±0.10 

3.71 

-19.0 

- 20.0 

-19.37 

14420 

0.0319 

0.1333 

4.45 ± 0.47 

1.76 ±0.09 

3.82 

-19.5 

- 20.5 

-19.82 

11122 

0.0397 

0.1500 

4.59 ±0.44 

1.76 ±0.07 

3.82 

- 20.0 

- 21.0 

-20.26 

4300 

0.0492 

0.1500 

5.52 ±0.88 

1.87 ±0.13 

3.46 

-20.5 

- 21.5 

-20.71 

1118 

0.0608 

0.1500 

6.33 ± 1.01 

2.01 ±0.29 

3.12 

- 21.0 

- 22.0 

-21.17 

198 

0.0749 

0.1500 

— 

— 

— 








